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The Euler-Maruyama scheme is known to diverge strongly and numerically weakly when applied to nonlinear 
stochastic differential equations (SDEs) with superlinearly growing and globally one-sided Lipschitz continuous 
drift coefficients. Classical Monte Carlo simulations do, however, not suffer from this divergence behavior of Euler's 
method because this divergence behavior happens on rare events. Indeed, for such nonlinear SDEs the classical 
Monte Carlo Euler method has been shown to converge by exploiting that the Euler approximations diverge only 
on events whose probabilities decay to zero very rapidly. Significantly more efficient than the classical Monte 
Carlo Euler method is the recently introduced multilevel Monte Carlo Euler method. The main observation of 
this article is that this multilevel Monte Carlo Euler method does - in contrast to classical Monte Carlo methods 
- not converge in general. More precisely, we establish divergence of the multilevel Monte Carlo Euler method 
for a family of SDEs with superlinearly growing and globally one-sided Lipschitz continuous drift coefficients. In 
particular, the multilevel Monte Carlo Euler method diverges for these nonlinear SDEs on an event that is not 
at all rare but has probability one. As a consequence for applications, we recommend not to use the multilevel 
Monte Carlo Euler method for SDEs with superlinearly growing nonlinearities. Instead we propose to combine 
the multilevel Monte Carlo method with a slightly modified Euler method. More precisely, we show that the 
multilevel Monte Carlo method combined with a tamed Euler method converges for nonlinear SDEs with globally 
one-sided Lipschitz continuous drift coefficients and preserves its strikingly higher order convergence rate from 
the Lipschitz case. 

1 Introduction 

We consider the following setting in this introductory section. Let T G (0, oo), d,m G N := {1, 2, . . . }, let (fi, T,V) be 
a probability space vifith a normal filtration {J^t)te[o.T] and let ^: $7 — M'^ be an J"o/'B(M'')-measurable mapping with 
E[||^||grf] < CXI for all p E [l,oo). Moreover, let /ii M"^ — > M'' be a smooth globally one-sided Lipschitz continuous 
function with at most polynomially growing derivatives and let a: M.'^ ^ jjdxm -^^ ^ smooth globally Lipschitz 
continuous function with at most polynomially growing derivatives. In particular, we assume that there exists a 
real number c S (0, oo) such that {x — y, ^(x) — ^J.{y))^d < c\\x — y||^d and \\<j{x) — a{y)\\^dxm < c\\x — yW^d for all 
x,y E M"^. These assumptions ensure the existence of an up to indistinguishability unique adapted stochastic process 
X : [0,T] X n ^ M.'^ with continuous sample paths solving the stochastic differential equation (SDE) 



for t G [0,T] (see, e.g., Alyushina pj. Theorem 1 in Krylov [27 or Theorem 2.4.1 in Mao [29]). The drift coefficient 
H is the infinitesimal mean of the process X and the diffusion coefficient a is the infinitesimal standard deviation of 
the process X. Our goal in this introductory section is then to efficiently compute the deterministic real number 
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where / : M*^ — > M is a smooth function with at most polynomiahy growing derivatives. Note that this question is 
not treated in the standard hterature in computational stochastics (see, for instance, Kloeden and Platen [26 and 
Milstein [30]) which concentrates on SDEs with globally Lipschitz continuous coefficients rather than the SDE ([ij. 
The computation of statistical quantities of the form ^ for SDEs with non-globally Lipschitz continuous coefficients 
is a major issue in financial engineering, in particular, in option pricing. For details the reader is refereed to the 
monographs Lewis |28j. Glassermann lOi, Higham [16 and Szpruch [37j . 

In order to simulate the quantity ([2]) on a computer, one has to discretize both the solution process X : [0, T] x £7 — > 
M'^ of the SDE ([T]) as well as the underlying probability space (r2,J^, P). The simplest method for discrctizing 
the SDE ^ is the Euler method (a.k.a. Euler-Maruyama method). More formally, the Euler approximations 
Y^-.n^ M'', 71 e {0, 1, ... , N}, N eN, for the SDE ^ are defined recursively through Fj^ := ^ and 

Yn+i Yn + • § + '^{Y^') - w^) (3) 

for all n G {0, 1, ... , A^} and all iV G N. Convergence of Eulcr's method both in the strong as well as in the 
numerically weak sense is well-known in case of globally Lipschitz continuous coefficients /i and a of the SDE (see, 
e.g.. Section 14.1 in Kloeden and Platen [35] and Section 12 in Milstein [30] )■ The case of superlinearly growing and 
hence non-globally Lipschitz continuous coefficients of the SDE is more subtle. Indeed, Theorem 2.1 in the recent 
article |22| shows in the presence of noise that Euler's method diverges to infinity both in the strong and numerically 



weak sense if the coefficients of the SDE grow superlinearly (see Theorem 2.1 below for a generalization hereof). 
In this situation. Theorem 2.1 in [22] also proves the existence of events VLj^ E J-, N G N, and of real numbers 
61, c e (l,oo) such that P^n] > 6'(-^') and |r/(w)| > c^'^") for aU uj e N e N. Clearly, this implies the 
divergence of absolute moments of the Euler approximation, i.e., limAr_j.co £[1^^?^ |^] — oo for aU p e (0, oo). 

The classical method for discrctizing expectations is the Monte Carlo Euler method. Let Y^^'': — >■ K"*, n G 
{0, 1, . . . , N}, iV g N, for fc e N be independent copies of the Euler approximations (|3| (see Section |3] for the precise 
definition). The Monte Carlo Euler approximation of ^ with iV e N time steps and Monte Carlo runs (see 
Duffie and Glynn [6 for more details on this choice) is then the random real number 

(4) 




Convergence of the Monte Carlo Euler approximations Q is well-known in case of globally Lipschitz continuous 
coefficients /i and a (see, e.g.. Section 14.1 in Kloeden and Platen [5B] and Section 12 in Milstein [30 ). Recently, 
convergence of the Monte Carlo Euler approximations Q has also been established for the SDE ([I]). More formally. 
Theorem 2.1 in [20] implies 
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P-almost surely (see also Theorem 3.1 below). The Monte Carlo Euler method is thus strongly consistent (see, e.g., 
Nikulin [33], Cramer [2] or Appendix A.l in Glassermann [TO])- The reason why convergence ^ of the Monte Carlo 
Euler method does hold although the Euler approximations diverge is as follows. The events J^at, TV S N, on which 
Euler's method diverges (see Theorem |2.1| below) are rare events and their probabilities decay to zero faster than 
any polynomial in as — > oo, sec Lemma 2.6 in [21] for details. Therefore for large N £ N the event fl^ is too 
unlikely to occur in any of Monte Carlo simulations in Q. 

Considerably more efficient than the Monte Carlo Euler method is the so-called multilevel Monte Carlo Euler 
method in Giles [5] (see also Creutzig, Dereich, Miiller-Gronbach and Ritter [3], Dereich [i], Giles [7], Giles, Higham 
and Mao [S], Heinrich [T31|T3], Heinrich and Sindambiwe [T? and Kebaier [531 for related results). In this method, 
time is discretized through the Euler method and expectations are approximated by the multilevel Monte Carlo 
method. More formally, let Y^^^''' : fl ^ W^, n G {0,1, . . . , N}, N e N, for / e Nq := {0, 1, 2, . . . } and /c G N be 
independent copies of the Euler approximations Js]) (see Section [6] for the precise definition) . Then the multilevel 
Monte Carlo Euler approximations for the SDE (HT are defined as 



fc=l 1=1 \k=l j 

for N G {2^,2^,2'^, . . .}. In the case of globally Lipschitz continuous coefficients of the SDE ([T]), this method has 
been shown to converge significantly faster to the target quantity ([2]) than the Monte Carlo Euler method Q. 
More precisely, in the case of globally Lipschitz continuous coefficients /i and u, the multilevel Monte Carlo Euler 
method ([6]) converges with order while the Monte Carlo Euler method converges with order |— with respect to 
the computational effort (see Section 1 in Giles [8] or Creutzig, Dereich, Miiller-Gronbach and Ritter 3 for details). 
In the general setting of the SDE ([T]) where /i does not need to be globally Lipschitz continuous, convergence of the 
multilevel Monte Carlo Euler method Q remained an open question. 
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Figure 1: Four sample paths of the approximation error of the muhilevel Monte Carlo Euler approximation Q for 
the SDE for iV e {2^, 2^, . . . , 2^} with T = 1. 



The convergence ^ of the Monte Carlo Eulcr method and the fact that Eulcr's method diverges on very rare 
events only shaped our first guess that the multilevel Monte Carlo Euler method should converge too. Surprisingly, 
however, convergence of the multilevel Monte Carlo Euler method fails to hold in the general setting of the SDE ([T]). 
To prove this it suffices to establish non-convergence for one counterexample which we choose to be as follows. Let 



m 



1, let fi{x) 



<j{x) = 0, f{x) = for all a; e M and let ^ : 17 — R be standard normally distributed. 



Clearly, this choice satisfies the assumptions of the SDE ([T]) and the SDE ([T]) thus reduces to the random ordinary 
differential equation 

dXt = dt, Xo = ^ (7) 

for t G [0,T]. The main observation of this article is that the approximation error of the multilevel Monte Carlo 
Euler method for the SDE ([t]) diverges to infinity. More formally. Theorem 
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below implies 
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P-almost surely. Note that the multilevel Monte Carlo Euler method diverges on an event that is not rare but has 
probability one. Thus - in contrast to classical Monte Carlo simulations - the multilevel Monte Carlo Euler method 
is very sensitive to the rare events on which Euler's method diverges in the sense of Theorem |2 . 1 1 below. To visualize 
the divergence ([s]) , Figure [l] depicts four random sample paths of the approximation error of the multilevel Monte 
Carlo Euler method ^ for the SDE Q with T = 1 and shows explosion even for small values of N & {2^, 2^, 2'^, . . .}. 
We emphasize that we are only able to establish the divergence ([s]) for the simple SDE ([t]) . Even in this simple case, 
the proof of the divergence ([S]) is rather involved and requires precise estimates on the speed of divergence of Euler's 
method for the random ordinary differential equation ([t]) on an appropriate event of instability; see below for an 
outline. 

Comparing the convergence result ^ for the Monte Carlo Euler method and the divergence result Q for the 
multilevel Monte Carlo Euler method reveals a remarkable difference between the classical Monte Carlo Eulcr method 
and the new multilevel Monte Carlo Euler method. The classical Monte Carlo Euler method applies both to SDEs 
with globally Lipschitz continuous coefficients and to SDEs with possibly superlinearly growing coefficients such as our 
SDE ([ij. The multilevel Monte Carlo Euler method, however, produces often completely wrong values in the case of 
SDEs with superlinearly growing nonlinearities. This is particularly unfortunate as SDEs with superlinearly growing 
nonlinearities are very important in applications (see, e.g., [28l|40l[37] for applications in financial engineering). We 
recommend not to use the multilevel Monte Carlo Euler method for applications with such nonlinear SDEs. 

Nonetheless, the multilevel Monte Carlo method can be used for SDEs with non- globally Lipschitz continuous 
coefficients when being combined with a strongly convergent numerical approximation method. Hutzcnthaler, Jentzen 
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and Kloeden [21] , for example, proposed the following slight modification of the Euler method ([s]) . Let 
n € {0, 1, ... , TV}, TV € N, be defined recursively through :— ^ and 
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for all n € {0, 1, . . . , A'' — 1} and all TV S N. Following [21] we refer to this numerical approximation as a tamed 
Euler method. Additionally, let Z^'^'^ : -> M, n e {0, 1, . . . ,TV}, iV G N, for ? e No and A: e N be independent 
copies of the tamed Euler approximations ([9]). In Theorem 6.2 below we then prove convergence of the multilevel 
Monte Carlo tamed Euler method for all locally Lipschitz continuous test functions on the path space whose local 
Lipschitz constants grow at most polynomially. In particular. Theorem |6.2| below implies the existence of finite 
random variables : — ?> [0, oo), e £ (0, |), such that 



E 



fix J 



1 

N 



N 



k=l 




< 



(10) 



for all N S {2^, 2^, 2^, . . .} and all e S (0, ^) P-almost surely. To sum it up, the classical Monte Carlo Euler method 
converges (see ([5])), the new multilevel Monte Carlo Euler method, in general, fails to converge (see (|8|) and the 
new multilevel Monte Carlo tamed Euler method converges and preserves its striking higher convergence order from 



the Lipschitz case (see (10)). Thus, concerning applications, the message of this article is that the multilevel Monte 
Carlo Euler method (|6]) needs to be modified appropriately when being applied to SDEs with superlinearly growing 
nonlinearities. This is a crucial difference to the classical Monte Euler method which has been shown to converge 
for such SDEs and which does not need to be modified. However, when modified appropriately (see, e.g., (|9|), the 
multilevel Monte Carlo method preserves its strikingly higher convergence order from the global Lipschitz case and 
is significantly more efiicient than the classical Monte Carlo Euler method even for such nonlinear SDEs. Thereby, 
this article motivates future research in the construction and the analysis of "appropriately modified" numerical 
approximation methods. 

For the interested reader, we now outline the central ideas in the proof of (Isl). For this we use the random 
variables : 
M, fc € N, I e 



n ^R, I eNa, k eN, defined by := Y^'''^''' for all M e N, / G No, /c G N. Then we note for every 



and every co & ft that |y^^''''^(a;)| is strictly increasing in n G {0, 1, ... , M} if and only if [^'^'^(w)] = 
It turns out that \Y^ {u!)\ increases in n € {0, 1, . . . , Mj d ouble exponentially fast for 



and Corollary 4.7 below for details). A central 



''''=(w)| > (2M)jr- 
aWuj G {l^'^'^l > (2Af)3r-3}, I G No and aU k,M G N (see Lemma[4.4 

observation in our proof of the divergence Q is then that the behavior of the multilevel Monte Carlo Euler method is 
dominated by the highest level that produces such double exponentially fast increasing trajectories. More precisely, 
a key step in our proof of (|8| is to introduce the random variables Ljy : — {1, 2, . . . , Id(A^)}, N G {2^, 2^, 2^, . . .}, 

by 

Ln :=max({l}u{;G{l,2,...,ld(7V)}: 3 fc G {1, 2, . . . , f } : > 2ir-i }) (11) 

for all N G {2\ 2^, 2^, . . .}. Using the random variables Ln, N G {2^2^,2^,...}, we now rewrite the multilevel 
Monte Carlo Euler method in (Isl) as 



^E(n''°ivT|(i:(v?'")'-(>3:':;;-' 



fc=i 



1=1 



, fc=i 



E 



/e{o,i,---,iog2(7V)} fc-i 



y 

iG{l,2,...,log2(JV)} k=l 



N 
2^N 



fc=l 



k=l 



2^«,Ljv,fc 



N 



2^N 

E/^2<^"-i),L„,fc 



(12) 
(13) 

(14) 



fe=i 



for aU N G {2^22,2^,...}. Due to the definition of Ljv, N G {2\ 2^, 2^, . . .}, it turns out that the asymptotic 



behavior of the multilevel Monte Carlo Euler method (12) is essentially determined by the three summands in (14) 



(see inequality (84|, estimate (91) and inequalities (96), (97) in the proof of Theorem 4.1 for details). In order 



to investigate these three summands, we - roughly speaking - quantify the value of the largest summand in each 



of the three sums in (14). For this we introduce the random variables rj^: [0, oo) and 9^:0,^ [0, oo) for 

iVG {21,22,23,...} by 

m :=max{|^^"'^-| gM: fcG {1,2,...,^}} (15) 

and 

9n :-max{|f^"-i^*^| gM: fcG {1,2,...,^^^}} (16) 
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for all TVs {2^ ,2'^ ,2^ , . . .}. Using the random variables 77 at : [0, 00) and 61^: [0,oo) for e {2^ ,2'^ ,2^ , . . .} 
we then distinguish between three different cases (see inequality ( p4| , inequality ( [9l] ) and inequalities (96), ( [QT] ) 

below). First, on the events {j]n > r~j} G J", iV g {2^,2^,2'^, . . . }, the middle summand in (14| will be 

positive with large absolute value and will essentially determine the behavior of the multilevel Monte Carlo Euler 

approximations ([l2]) (see estimate (84) for details). Second, on the events {rjN < 2 'i T~i} D {rjN < 6^} e J^, 
N £ {2^, 2^, 2'^, . . . }, the left summand in ( |l4| ) will be positive with large absolute value and will essentially determine 
the behavior of the multilevel Monte Carlo Euler approximations (12) (see inequality (91 1 for details). Finally, on 

the events < 2^^t^T"3} n {77Ar > On} e J^, N e {2\ 2^, 2^, . . . }, the right summand in (fl4| wih be 

negative with large absolute value and will essentially determine the behavior of the multilevel Monte Carlo Euler 
approximations (12) (see inequalities (96) and ( [97| for details). This very rough outline of the case- by-case analysis 
in our proof of ([8|) also illustrates that the multilevel Monte Carlo Euler approximations ( 12 ) assume both positive 
(first and second case) as well as negative values (third case) with large absolute values. We add that this case-by-case 
analysis argument in our proof of (jsj) requires that the probability that the random variables rj^ and 9^ are close 
to each other in some sense must decay rapidly to zero as g {2^, 2^, 2'^, ... } goes to infinity (see inequality ( 71 ) 
below). We verify the above decaying of the probabilities in Lemma 4.13 below which is a crucial step in our proof 
of ([S]). Additionally, we add that the level Ljv is approximately of order log(log(iV)) as N goes to infinity (see 
Lemma 4.9 for the precise assertion). In view of the above case- by-case analysis of the multilevel Monte Carlo Euler 
method, we find it quite remarkable to observe that the essential behaviour of the multilevel Monte Carlo Euler 
method in (|8| is determined by the levels around the order log(log(A^)) as N goes to infinity. 

The remainder of this article is organized as follows. Theorem |2.1| in Section [2] slightly generalizes the result on 
strong and weak divergence of the Euler method of Hutzenthaler, Jentzen and Kloeden [22 . Convergence of the 
Monte Carlo Euler method is reviewed in Section |3] The main result of this article, i.e., divergence of the multilevel 
Monte Carlo Euler method for the SDE ([7|, is presented and proved in Section[4] We believe that the multilevel Monte 
Carlo Euler method diverges more generally and formulate this as Conjecture |5.1| in Section [5j Section |6] contains 
our proof of almost sure and strong convergence of the multilevel Monte Carlo tamed Euler method for all locally 
Lipschitz continuous test functions on the path space whose local Lipschitz constants grow at most polynomially. 



2 Divergence of the Euler method 

Throughout this section assume that the following setting is fulfilled. Let T E (0, 00), let {n,J-,P) be a probability 
space with a filtration {J^t)telo,T] and let : [0, T] x f2 — >■ K be a one-dimcnsional standard (J-'f )tg[o,T]-Brownian 
motion. Additionally, let ^: 17 — >■ M be an Jx)/'B(M)-measurable mapping and let cr : M — >■ M be two B{M.)/B{M.)- 
measurable mappings. We then define the Euler approximations : 11 — > M, n e {0, 1, . . . , N}, N £N, recursively 
by := e and 

Yn+i ■■= + KYn) ■ § + ^Y^:') ■ - W.^) (17) 

for all n € {0, 1, . . . , — 1} and all iV e N. The following theorem generalizes Theorem 2.1 in Hutzenthaler, Jentzen 
and Kloeden [22] . 

Theorem 2.1 (Strong and weak divergence of the Euler method). Assume that the above setting is fulfilled and let 
a,c Cz (1,00) be real numbers such that \fJ-{x)\ + |cr(a;)| > -1^^ for all x € M. with \x\ > c. Moreover, assume that 
P[cr(^) 7^ 0] > or that there exists a real number (3 G (l,oo) such thatV\\£\ > x] > I3^~^^^ for all x G [l,oo). Then 
there exists a real number 9 £ (1, 00) and a sequence of nonempty events 17 jy (z J-, N CI N. such that P[l7jv] ^ 9^^^ ) 



and, for all UJ e 17 jv and all N £ N. In particular, the Euler approximations ( |17[ ) satisfy 

limA,^ooIE[|r^|^'] = 00 for allpe {0,oo). 



Theorem |2 . 1 1 immediately follows from Lemma \2?2\ and Lemma [273| below. More results on Euler's method for SDEs 
with possibly superlinearly growing nonlinearities can, e.g., be found in [12U111I5T1I32] and in the references therein. 

Lemma 2.2 (Tails of , N G N). Assume that the above setting is fulfilled and let P[it(^) 7^ O] > 0. Then there 
exists a real number (3 G (l,oo) such that V[\Y-^\ > x] > /3(-(^^)'') for all X G [1, 00) and all iV G N. 



Proof of Lemma 2.2 By assumption we have P[|(t(^)| > O] > 0. Therefore, there exists a real number K G (l,oo) 
such that 



g(0,cx)). (18) 
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Moreover, we have 



> 



> 



|C + KOw + <^iOWT I > x\ >v[\aiOWT I - lel - T|/i(e)| > X 

Wm > ^, 1^1 + THOl < K, \a{OWT I - lei - r|M(OI > X 

> j?, m + T\KO\ <K,j^\WT\-K> 



for all X € [l,c«) and all N ^N. Definition (18) and Lemma 4.1 in (35] therefore show 





> p 







K ' 



IS- \ \ 



T"2 \Wt\ > T-^N^Kix + K) 



> — = • expf - T-^NK^(x + K)^) > — = • expf - AT-^K'^iNxY 



4VT 



> e 



4VT 



for all X e [1, oo) and all N £N. This completes the proof of Lemma [2^ 



□ 



Lemma 2.3. j45SMme i/iat the above setting is fulfilled and let a, c e (1, oo) be real numbers such that \fi{x)\ + \a{x)\ > 
for all X Cz M. with \x\ > c. Moreover, assume that there exist real numbers Nq G {0, 1,2,...}, /3 e (1, oo) such 
that V[\Y^J >x]> for all x £ [1, oo) and all N eNn {Nq, iVo + 1, . . . }. Then there exists a real number 

9 S (1, oo) and a sequence of nonempty events G J^, G N H {Nq, No + 1, . . . }, such that P[i7jv] > 9^^^ ) and 



Y^ico) >c 



for all oj E and all N £ N {Nq, Nq + 1, . . . }. In particular, the Euler approximations 



sails/?/ limAr_j.oo E [I Y^|P] — oo for all p £ (0, oo). 

Proof of Lemma\2.3\ Define real numbers rjv G [0, oo), G N, by 



tn '■= max c. 



(19) 



for all G N. We also use the function sgn: M — > M defined by sgn(x) := 1 for all x G [0, oo) and by sgn(a;) := —1 
for all X G (— oo, 0). Furthermore, we define events G J", iV G N n {Nq,Nq + 1, . . . }, by 



■■= (^at-^ jc. G 17 : sgn(M(i;f (c.)) • a{Y^ {lo))) ■ {Wi^{u:) - W^{u:)) > 



flLGf7:|r/„H|>(r^) 



((^) 



(20) 



for all G N n {Nq, A^o + 1, • ■ • }■ In particular, the definition of {^N)NeNn{No.No+i,... } implies 



T 
N 



/i(r„^^(^)) + (u;)) • Wi^iuj) - W^i^) 



for all n G {A^o, A^o + 1, . . . , Af - 1}, a; G and aU A^ G N n {A^o, A^o + 1, • • ■ }• 

In the next step let A^ G N n {Ao, A'o + l, ■ • • } and u E he arbitrary. We then claim 

|r,r(c.)|>(r^)«"*')") 



(21) 



(22) 



for all n G {Nq, Nq + I, . . . ,N}. We now show ( 22 ) by induction on n G {Nq, Nq + 1, . . . , N}. The base case n = Nq 
follows from definition (20 1 of ri^v- For the induction step assume that ( 22 ) holds for one n G {Nq, Ao + 1,...,A^ — 1}. 
In particular, this implies 



K{u^)\ > (rA.)«^)") 



> rpf > c > I. 



(23) 



Moreover, definition (17), the triangle inequality and equation (21) yield 

T 



T 

N 



T 



^(y,f (c.))| + |a(r„^(c.))| . \w^{u)-W^{u;)\ - (a.)| > -(|^(y,f (a.))| + \a{Y„^{u))\)- |y„^(a.) 
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and the estimate |/i(a::)| + > ^j— for all a; G M with |x| > c, inequality (23) and definition (19 1 therefore show 

, „ , (° + i) / T , „ , \ , „ , (°+i) / T \ , , (°+i) 

The induction hypothesis hence yields 



Inequality ^ thus holds for ah n e {A^o, iVo + 1, . . . , N}, uj e and all iV e Nn{iVo, iVo + 1, ■ ■ • }■ In particular, 
we obtain 

\Y^i^)\ > (r-iv)(^'^^") > ci^^^") (24) 
for all a; G r^AT and all iV G N n {Nq, Nq + I,. ■ ■}■ Additionally, Lemma 4.1 in [55] yields 



E 



^|sgn(M(y„")-^(F„"))-(^ty („H^i)T -w^)>i| 



(25) 







T 




)>- 


= P 


Wt > — 


= P 


/ - N_ 




N - N 





> 



N 



'-almost surely for all 7i G {0, 1, . . . , — 1} and all iV G N. Therefore, we obtain 

^ x(N-No) 



T—^Wt > V ^ 

AT 



>/3 



(-(JVr„)((°2') 



> e 



for all G N n {Aq, A^q + 1, . . . }. This shows the existence of a real number 9 G (1, oo) such that 
for all AT G N n {A^o, + 1, • ■ • }■ Combining ([24| and ([26| finally gives 



(26) 



hm E 

N—yoo 



\VN\P 



> lim E 



1 \V^\P 



> lim ( P^r^Ar] •cl^'^"^')") ) > Jim 



-TV 



□ 



for all p G (0, oo). This, (24) and (26) then complete the proof of Theorem 2.1 

3 Convergence of the Monte Carlo Euler method 

The Monte Carlo Euler method has been shown to converge with probability one for one-dimensional SDEs with 
globally one-sided Lipschitz continuous drift and globally Lipschitz continuous diffusion coefficients (see dO])- The 
Monte Carlo Euler method is thus strongly consistent (see, e.g., Nikulin [33], Cramer [5] or Appendix A.l in Glasser- 
mann '10]). After having review this convergence result of the Monte Carlo Euler method, we complement in this 
section this convergence result with the behavior of moments of the Monte Carlo Euler approximations for such 
SDEs. More precisely, an immediate consequence of Theorem |2.1| is Corollary |3 . 2 1 below which shows for such SDEs 
that the Monte Carlo Euler approximations diverge in the strong L^-sense for every p G [l,cx)). We emphasize that 
this strong divergence result does not refiect the behavior of the Monte Carlo Euler method in a simulation and it 



is presented for completeness only. Indeed, the events on which the Euler approximations diverge (see Theorem 2.1) 
are rare events and their probabilities decay to zero very rapidly (see, e.g.. Lemma 4.5 in [2CP for details). This is 
the reason why the Monte Carlo Euler method is strongly consistent and thus does converge according to [20] (see 
also Theorem 3.1 below). 

Throughout this section assume that the following setting is fulfilled. Let T G (0, oo), let (fi, J^, P) be a probability 
space with a normal filtration {J^t)telo,T], let W'' : [0, T] x — M, fc G N, be a family of independent one-dimensional 
standard (J^t)fg[g j^j-Brownian motions and let ^'^ : — > M, fc G N, be a family of independent identically distributed 
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J-o/S(K)-nieasurable mappings with E [1^1 |P] < oo for all p E [1, oo). Moreover, let /x, : M ^ M be two 
measurable mappings such that there exists a predictable stochastic process X : [0, T] x J7 -> M which satisfies 
/J" |/i(Xs)| + |cr(X3)|^ ds < oo P-almost surely and 



(27) 



P-almost surely for all t E [0,T]. The drift coefficient fi is the infinitesimal mean of the process X and the diffusion 
coefficient a is the infinitesimal standard deviation of the process X. We then define a family Y^''^ : — M, 
ri e {0, 1, . . . , N}, N,k E N, of Euler approximations by Fq ' :— ^'^ and 



N,k 
n+1 



(28) 



for all n G {0, 1, . . . , — 1} and all N,k E N. For clarity of exposition we recall the following convergence theorem 
from |5D]. Its proof can be found in [5D]. 

Theorem 3.1 (Strong consistency and convergence with probability one of the Monte Carlo Euler method). Assume 
that the above setting is fulfilled, let ^,a, f : M — M he Jour times continuously differentiable and let c E [0, oo) be a real 
number such that {x — y) ■ (/i(a;) — fi{y)) < c |x — y| , l^^ix) — c(y)| < c jx — j/| and \fi^'^'^{x)\ + |cr('*)(a;)| + [/'■^^(a;)| < 
c(l + \x\ ) for all X eR. Then there exist finite J- /B{[0, oo))-measurable mappings : il — > [0,oo), e E (0,1), such 
that 



E 



> fc=l 



< 



(29) 



for all N E N and all e E (0, 1) V-almost surely. 



In contrast to pathwise convergence of the Monte Carlo Euler method for SDEs with globally one-sided Lipschitz 



continuous drift and globally Lipschitz continuous diffusion coefficients (see Theorem 3.1 above for details), strong 



convergence of the Monte Carlo Euler method, in general, fails to hold for such SDEs which is established in the 



following corollary of Theorem 2T i.e., in Corollary |3.2[ As mentioned above we emphasize that Corollary 3^ does 
not reflect the behavior of the Monte Carlo Euler method in a practical simulation because the events on which the 



Euler approximations diverge (see Theorem 2.1) are rare events and their probabilities decay to zero very rapidly 
(see Lemma 4.5 in [20 for details). 

Corollary 3.2 (Strong divergence of the Monte Carlo Euler method). Assume that the above setting is fulfilled and 
let a,c E (l,oo) be real numbers such that \fJ-{x)\ + \(7{x)\ > ^j— for all x E M. with \x\ > c. Moreover, assume that 
P[a{£_^) 7^ O] > or that there exists a real number (5 E (l,oo) such that P[|^ for all X E [1, oo). 

Moreover, let f : R ^ R be B{R)/B{R) -measurable with f{x) > l\x\^ - c for all x E R. Then 



lim E 



E 



fiXr) 




(30) 



for all p E [1, oo). 



Proof of Corollary \3.^ The triangle inequality, Jensen's inequality and the estimate f(x) > ^ja^l^ c for all x E 
give 



E 



fiXr) 



7V2 



, fe=i 



> 



iV2 



Lp{n-M) 



N , 



fc=l 



E 



1/(^1 



LP(0;B 



> ^E 

- 



k = \ 



(31) 



E 



|/(^t)| 



= E 



!{Yi 



Af.l> 

N , 



-E 



1/(^1 



> - -E 

c 



Y, 



N 



E 



|/(^t)| 



for all G N and all p E [l,oo). Combining (31) and Theorem 2.1 then shows (30l in the case E[|/(X7-)|] < oo. In 
the case E[|/(X7-)|] = oo, the estimate f{x) > — c for all x G M shows E[/(X7-)] = oo and this implies (30) in the 
case E[|/(Xj-)|] = oo. The proof of Corollary 3.2 is thus completed. □ 
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4 Counterexamples to convergence of the multilevel Monte Carlo Euler 
method 

Theorem |4.1| below establishes divergence with probability one of the multilevel Monte Carlo Euler method Q 
for the SDE ([7|. This, in particular, proves that the multilevel Monte Carlo Euler method is in contrast to the 
classical Monte Carlo Euler method, in general, not consistent (see, e.g., Nikulin '37, Cramer or Appendix A.l 
in Glassermann [ID]). 

Throughout this section assume that the following setting is fulfilled. Let T,a G (0,00), let {^l,J^,¥) be a 
probabihty space and let ^'■'^ : — >■ K, Z G No, S N, be a family of independent normally distributed J^/B{R)- 
measurable mappings with mean zero and standard deviation a. Moreover, let X : [0, T] x 51 — > E be the unique 
stochastic process with continuous sample paths which fulfills the SDE 



dXt = -X^ dt, Xo^^ 

N,l,k . 



for t G [0,T]. We then define a family of Euler approximations y,j ' ' : 51 
fc e N, by Yr 



N.Lk 



and 



N,Lk 
n+1 



NJ.k 



V n 



N.l,k\^ 

' N 



(32) 

71 e {0,1,...,7V}, iV e N, Z e No, 

(33) 



for all 71 e {0, 1, . . . , iV - 1}, e N, / e No and all fc e N. 

Theorem 4.1 (Main result of this article: Divergence with probability one of the multilevel Monte Carlo Euler 
method for the SDE (32)). Assume that the above setting is fulfilled. Then 



lim 

ld(Af)eN 



fc=l 



A^ 

1=1 k=i 



1) I 



(34) 



-almost surely for all p e (0,cx)). 
The proof of Theorem |4.1| is postponed to Subsection |4.2| below. 



4.1 Simulations 



We illustrate Theorem 14 . II with numerical simulations. To this end we observe that the exact solution of the random 



ordinary differential equation ( 32 ) satisfies 



Xt = 



(35) 



J 

[l + At^^Y 

for all t e [0, 1]. The real number E[(Xi)2] can then be computed by numerical integration or by the Monte Carlo 
method. Figure [l] depicts four random sample paths of the approximation error of the multilevel Monte Carlo Euler 
approximations in the case T = 1 and ct = 1 in (32) where E[(Ai)2] 0.28801. The sample paths clearly diverge 
even for small A^ e {2^, 2^, 2^, . . .}. For some other SDEs, however, pathwise divergence does not emerge for small 
A^ e {2^, 2^, 2^^, . . .}. For example, let us choose a standard deviation as small as ct = 0.1 in (32) where T = 1. Here 
the exact value is E[(Ai)2] « 0.009971. Then sample paths of the multilevel Monte Carlo Euler approximation seem 
to converge even for reasonably large A^ g {2^, 2^, 2"^, . . .} (see Figure|2]for four sample paths). So the sample paths 
of the multilevel Monte Carlo Euler method for some SDEs first seem to converge but diverge as A^ G {2^, 2^, 2^^, . . .} 
becomes sufficiently large. To see this in a plot, we tried different values of a and found sample paths in case of 
a = I and T = 1 which first seem to convergence to the exact value E[(Xi)2] « 0.09248 but diverge for larger values 
of A^ G {2^2^, 2^, . . .} (see Figure [3]for four sample paths). 



4.2 Proof of Theorem liTTl 



First of all, we introduce more notation in order to prove Theorem 
X e M, be defined recursively through ''^ := x and 

.5 T 



4.1 



N,x _ N,'. 
Vn+l '■— Vn 



N.L 



1 



Let y^'- e 



.N' 



[Vn 



' N 



e {0,1,..., A^}, A^ e N, 



(36) 



for all n E {0, 1, ... ,A^ — 1}, A^ G N and all a; G K and let p G (0,oo) be fixed for the rest of this section. This 
notation enables us to rewrite the multilevel Monte Carlo Euler approximation in (34) as 



1 ^ 



1,0, fe IP 



k=l 



1=1 k=i ^ 



Y 



2"-'^\Lk\P 



2('-i) 



Id(Ar) 

1=0 k=l 



A^' 

1=1 k=i 



ld(W) ST 
1=0 k=l 



1=1 k=l 



(37) 
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10" 



CD 

o 



J5 10" 



X 

o 

a 10- 




10" 



10' 

N 



10° 



10"' 



10" 



10"' 




10" 



10' 

N 



10° 



10" 



g 

I 10" 

■>< 
o 

I 10" 
< 




10 



10 

N 



10 



10"' 



10"' 



10"' 




10" 



10 

N 



10 



Figure 2: Four sample paths of the approximation error of the multilevel Monte Carlo Euler approximation in (34) 
where T ^ 1, a = 0.1, p = 2 and N £ {2^,2'^, 222}. 



10= 



o 



« 10" 



X 

o 

a. 
a. 
< 



10" 




10" 



10' 

N 



10° 




10' 



o 

■« 10" 
i 
x 
o 



< 



10" 
10"' 



10" 



■ MLMC Euler 
Order line 1/2 



10' 

N 



10° 



10" 



10" 



10"' 



10" 




10' 

N 



10° 



Figure 3: Four sample paths of the approximation error of the multilevel Monte Carlo Euler approximation in (34) 
where T= 1, ct = i p = 2 and G {2^ ,2^ , . . . ,2^^}. 
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for all N G {2\ 2^, 2^, . . . }. Additionally, let Ln : ^ ^ {1,2, ... , Id(iV)} be defined as 

:=max({l}u{;e{l,2,...,ld(7V)}: 3 k e {1,2, ...,§} : jC'^'^l > 2ir-i } 
for every N £ {2^2^, 2^, . . .}. Furthermore, define 1]^ : ft ^ [0, oo) and 6n-^^ [0, oo) by 

:=niax{|^^"^'=| eM: fc e {l, 2, . . . , ^} } 



and 



iN max 



{|e(^«-i)''=|GM:fce{l,2,...,5^}} 



(38) 



(39) 



(40) 



for every N £ {2^,2^,2^, . . .}. Moreover, we define the mappings [•], [-J : M ^ Z by \x~\ := min{z E Z: z > x} and 
by \_x\ := max{z G Z: z < x} for all a; G M. Additionally, we fix a real number S G (0, |) for the rest of this section. 



In the next step the following events are used in our analysis of the multilevel Monte Carlo Euler method. Let A^^', 



A^^\ A^^\ A^^^ gT, N e{2\2\2^...}, be defined by 



:= < [2ld{<j'Thn{N))\ } (41) 

:={3;G{0,l,2,...,ld(7V)}3A:G{l,2,...,f}: l^''"] > 2^T-i n] (42) 
A^^= {3 ; G N, [2 \d{a^Thn{N))\ < I < ld{N) + 1 : 2JT-3 < ?7jv < 2JT-3 (l + 5(-'5-2<'-'')) | (43) 



A 



(4) 
TV 



(44) 



for aU N G {21,2^23,...}. Additionally, define Nq G {2^, 2^, 23, . . . } and Ni: fl ^ {2\ 2^, 2^, . . . } U {cx)} by 



Nn 



and by 



iVi(w) := min({c5o} U |n G {A^o, 2^Aro, 2^Ao, . . . } ; Vm G {n, 2^^, 2^71, . . .} : uj ^ A^^^ U A^^^ U A^^^ U A^^^j) (45) 
for all w G ri. Next we prove a few lemmas that we use in our proof of Theorem |4.1[ 

Lemma 4.2 (Dynamics for small initial values). Assume that the above setting is fulfilled. Then we have \yn'^\ ^ 
|a;| < (^) ^ for all n G {0, 1, . . . , A^} and all N eN. 

Proof of Lemma Fix A^ G N and |a;| < (^)^- We prove jy^'^^l < |a;| by induction on n G {0, 1, . . . , A^}. The 
base case n = is trivial. For the induction step n — )• n + 1, note that the induction hypothesis implies 

T 



N,x 
Vn+l 



.N,i 



Vn 



1 



N 



for all n G {0, 1, . . . , A^ — 1}. This completes the proof of Lemma 4.2 



(46) 



□ 



Lemma 4.3 (Dynamics for large initial values). Assume that the above setting is fulfilled. Then we have jj/^'^] > 
> (^)^ for all n G {0, 1, . . . , A^}, \x\ > (^)^ and all N eN. In particular, we have 



N,x 
Vn+l 



(47) 



for alln(^{Q,l,...,N - 1}, \x\ > (^) " and all A G N. 



Proof of Lemma\l^ Fix A^ G N and |x| > (^) ^ We prove |?/^'"^| > \x\ by induction on n G {0, 1, . . . , A^}. The 
base case n = is trivial. For the induction step n — >■ n + 1, note that the induction hypothesis implies 



N.x 

Vn+l 



N,x\ 



(48) 



for all n G {0, 1, . . . , A^ — 1}. This completes the induction. The assertion (47) then immediately follows by taking 
absolute values in (36 1. □ 



Lemma 4.4 (Growth bound for large initial values). Assume that the above setting is fulfilled. Then we have 

(5") 



'T\4 I AT.^I 



< 



((^)^n) 



(49) 



for all n G {0, 1, . . . , N), \x\ > (^) " and all N e N 



11 



Proof of Lemma\l4\ Fix G N and |a;| > (^) ^ . We prove ( |49[ ) by induction on rt G {0, 1, . . . , N}. The base case 
n = is trivial. For the induction step n — >■ n + 1, note that Lemma |4.3| and the induction hypothesis imply 



N.x 
Vn+l 



< 



(5") 



for all n G {0, 1, . . . , — 1}. This completes the proof of Lemma 4.4 



Lemma 4.5 (Monotonicity) . Assume that the above setting is fulfilled. Then we have 



> 



(50) 

□ 

(51) 



for alln E {0,1,...,^}, all x,y E R satisfying \x\ > \y\, \x\ > i^)" and all N E N. 



Proof of Lemma \4-5\ Fix N E N and x,y E R with \x\ > \y\, \x\ > {^)'^- We prove (51 1 by induction on 
n E {0,1,..., N}. The base case n = is trivial. For the induction step n — > n + 1, note that Lemma 4.3 and the 
induction hypothesis imply 



Vnil 



.,N,x\ 



T 
N 



1 > b: 



N,y\ 



T 

N 



.,N,x\ 



N,v\ 



T 
N 



.,N,y\ 



1 



N,y 
Vn+l 



for all n E {0, 1, . . . , iV — 1}. This completes the proof of Lemma 4.5 



(52) 



□ 



Lemma 4.6 (Dynamics of multiples of the initial value). Assume that the above setting is fulfilled. Then we have 



for allnE {0, 1, . . . , TV}, > M E [l,oo) and all N eN. 



(53) 



Proof of Lemma \4^ Fix N E N. We prove (53) by induction on n G {0, 1, ... , N}. The base case n = is trivial. 
For the induction step n — > ri + 1, note that Lemma |4.3| and the induction hypothesis imply 



N,Mx 
Vn+l 



> 



- 1 



N,x 
Vn+l 



for all n E {0, 1, . . . , iV — 1}, |x| > (4^) * and all M E [1, oo). This completes the proof of Lemma 



4.6 



(54) 



□ 



Corollary 4.7. Assume that the above setting is fulfilled. Then we have \yn'^\ ^ )(^)* for all n E 

{0,1,.. .,iV}, |a;| > M(^)^ M E [l,oo) and all N E N. 



Proof of Corollary \4.T\ Lemma 



4.5 



Lemma 



4.6 



> 



Vn 



and |y, 



(4^)i imply 



y-n 



for all n G {0, 1, . . . , iV}, |a;| > M{^y , M E [1, oo) and all N eN. This completes the proof of CoroUary 
Lemma 4.8. Assume that the above setting is fulfilled. Then we have 

1 fr(l-r)N\ 



4.7 



N,x 

Vn 



(55) 



□ 



(56) 



for all \x\ > (^)' (1 + 5(-''^)), TV G N and all r E (0,oo). 

Proof of Lemma \4.8\ We apply the inequality 1 + z > exp(|) for all z E [0, 2]. Noting that 5^^''^-' < 1 < 2 for all 



iV G N and all r G (0, oo), we infer from Corollary 4.7 



N,2 

Vn 



(^)^(l + 5(-^-')*' ^(^)^[exp(i5(--))]^' ' = 
for all \x\ > (^) ^ (1 + 5("''^)) , iV G N and aU r E (0, oo). This completes the proof of Lemma 



4.8 



(57) 



□ 



12 



Lemma 4.9. Assume that the above setting is fulfilled. Then we have 



oo 



A. 



(1) 



< oo. 



(58) 



Proof of Lemma\4^ The definition ([38]) of L^, N e {2\ 2^, 2^, . . . }, and independence of I e Nq, k e imply 
pU^^j ^ P[V/ e N, [2ld{a^Thn{N))\ < I < ld(iV)Vfc e {1, 2, . . . , f} : < 2^T- 



Id(JV) 



ld(7V) 

n 



Id(Af) 

n 



|^0,1| < 24^- 



1 



;=L2 1d(ff2\/Tln(JV))J 

for all N e {A^o, S^A^o, 2^A^o, ■ ■ ■}• The inequality 



1 



> 



27r 



a;^ / - - a; I e 2.2 



27r 

3- V2)e-J" 



e^^dy > 2 



5 1 



a;e <i 



TT (73 + 72) ~ 6 



e 2 (fy 



> -xe ^ 



for all X € [0, 00) therefore yields 



ld(N) / 

^ n 

/=[2 1d(ff2VTln(JV))J V 



1 2i /'_3 24 ' 



< 



< 



< 



< 



1 



2L21d(ff2yTln(Af))J/4 



1 - 



1 



6aTt 
1 



• exp 



exp 



). « 
\ 2L2 ld(&2^,„(„))J 



3 . 2l-21d(5-^yTln(Ar))J/2\ \ 2 [2 id(<T2 ^i„(jv))j 



1 ^ • exp 

1 ^ • exp ( 

l-^-exp(-XA^) 



3 . 2id{S^VTln(N)) \ \ 2L2id(<72vTi„(N))J 



4a VT 



3-a-2Vrin(A^) \ \ 2 



2 1d(CT^\/Tln(iV)) 



4ct2Vt 



2ld((a2^1„(jV))2) 



= 1- 



1 



6aT i 



■ N- 



(0-2 vnn(lV))2 



for all iV e {iVo, 2iA^o, S^A^o, • ■ ■}• Next we estimate 1- x < exp(-x) for all x e M to get 



(59) 



(60) 



OO 



A. 



(1) 



E 



A'e{2^2^,2^,...} 
N<No 



A 



E 



Ne{2^ ,2^,2^,...} 
N>No 



A 



<No + 



E 



JVe{2^2^,2^,...} 

A'>JVo 



1 3 \ (52^1„(W))2 

1 ^-N-- 



<iVo+ E 



N=No 



exp 



(ff2^1„(jV))2 



iVo+ E 



N=Na 



6CT5Ti(ln(7V))2 



< oo. 



This completes the proof of Lemma 4.9 



Lemma 4.10. Assume that the above setting is fulfilled. Then we have 



C30 

E^I 



A. 



(2) 



< OO. 



□ 



(61) 
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Proof of Lemma \4-10\ Subadditivity of the probability measure P and the inequahty P[f7 ^|^'''"'^| > x] < ^ exp (— ^) 
for aU X € (0,00) (e.g., Lemma 22.2 in [24 ) imply 



A 



(2) 



N 
ld{N) 



3Ze {0,i,2,...,id(iV)}3fce{i,2,...,f}: l^'^'^l > 2^r'iiv' 



Id(JV) 



< 



1=0 k=l 



ld{N) 



1=0 



1=0 



2' 



a-'\e''\> 



< -^exp 



2ct2T5 



(Id(iV) + 1)ct2iT3 exp - 



(i-i) 



\d(N) _ 1 

7;r ■ , n-i> exp 



C-i) 1 

> 2 — T-iN 



(i-i) „^ 



2' 2^iV 



2a^Ti 



2-2(j^T-2 



for ah TV e {2\ 2^, 2^, . . .}. Summing over N e {2\ 2^, 2^, . . .} resuhs h 



^ p[a^^„^ < (Id(iV) + 1) CT2JT3 exp I ^— r I < 00 

n=l Are{2i,22,23,...} 

and this completes the proof of Lemma [4.10| 

Lemma 4.11. Let Z : CI ^ R be a standard normally distributed J- / B{M.) -measurable mapping. Then 



\Z\<x + y \Z\>x 



< bxy 



(62) 

(63) 

□ 

(64) 



for all X ^ [| , 00) and all y G [0, 00) . 

Proof of Lemma \4-ll\ Monotonicity of the exponential function yields 



[x <\Z\ < X + y]^2-V[x < Z < x + y] ^2 



x+y 



1 



e ~2 dz < 



ye 



(65) 



V27r V27T 

for all x,y [0, 00). Apply the standard estimate V[\Z\ > x] > -^^^2 exp ( — ^) for all x e (0, 00) (e.g., Lemma 



22.2 in 24J), inequality ([65]) and > ^ for all a; e [^,00) to get 



\Z\<x + y \Z\>x 



[x< \Z\ <x + y] 
V[\Z\ > x] 



< 



.ye 



xy 



< 5xy 



for all X e 00) and all y £ [0, 00). This completes the proof of Lemma 
Lemma 4.12. Assume that the above setting is fulfilled. Then we have 



4.11 



00 



A. 



(3) 



< 00. 



Proof of Lemma\4J^ Let : Nq x No x O ^ N U {od} be defined as 

K{v,l) := min (|fc G N: > 2iT-i| U {00}) 

for all v,l £ Nq. Inserting definition (|43| we get 



□ 



(66) 



(67) 



l(3) 



< 



< 



< 



'[3/ G N, [2ld{a^Thii{N))\ < I < Id(iV) + 1; 2iT-J <m < 2Jr-i (l + S'^^'^'' ")) 
-1 

pMr-i <r,w <2ir-i(i + 5(-*-2""'))" 



ld(Af) + l 

E 

l=[2 ld(ff2\/Tln(Af))J 
ld(Af) + l 

y phi; e {l,2,...,ld(7V)}: 2ir-i < max < 2iT~i (l + S'-^''^*'" 

;=[21d(&2v/Tln(Af))J ^ 
ld(W) + l ld(7V) 

^ ^ p[{3fcG{i,...,f}: |r'1>2^r-i}n 



(68) 



;=[21d((j2v/Tln(Af))J ^=1 

{vfce{i,...,f}: <2Jr-i(i + 5(-*-2'"'') 

ld(Af) + l Id(Af) 

< E Eip[{^(-'0<i^}n{|e 

;=[21d(CT2\/Tln(Af))J "=1 



ti,if(i;,;) 



< 2iT"i(l + 5 



(-A--2('-i)) 
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for all N G {iVo, 2^A'o, ^^No, . . .}. The method of rejection sampling results in 



A 



ld(Af) + l Id(Af) 
< ^ ^P[|r 

l=l2\d{a^VT\n{N))\ 

ld{N} + l ld(N) 

Z=[21d(5-2yTln(Af))J i'=l 
ld{N) + l 



K{v,l) 



),K{v,l) < oo 



< 



2ir-i(i + 5(-*-2""'')) 



\d{N) ■ P 

l=[2 ld(a^VTln(N))\ 



24 1+5 



(-^.2(^-l))^ 



1^0,11 > 2iT- 



-ltO,l| ^ 24 



5T4 



for all TV e {iVo, 2^A'^o, 2^A'^3, . . .}. In order to apply Lemma 



4.11 



we note that 



2^ 2L2'^(^^^'"(^))J/'* 2(2''^(*^^'"(^))^i)/'* 2'd(^^^'"(^))/2 
> , > , = 



crTi 



(tT i 



o-ri23 



> — ^ > 
2i - 21 " 2 



(69) 



for ain e N n [[21d(a-2T3ln(7V))J , oo) and aU N e {2^2^, 2^, . . .}. Lemma 4.11 applied to the standard normally 
distributed variable tT~^^°'^ thus leads to 



l(3) 



< ld{N) 



ld(Af) + l ld(Af) + l 

y 5 . ^ . ^5(-5-2<'-') ^ id(iv) y _ ^ 

/= [2 ld(5-2yT ln(Ar))J /= [2 ld[a^VTln{N))\ 

< Id(iV) ''y^' ^2i^^(_..2(.M(.^^M«,-i,) ^ (ld(iV))2i^5(-|-('^^^'"W)^) (70) 

/=[2 ld(5-2yTln(Af))J 

<l^expf-^^^^^:M^^ 
" a^T2 M 4 i 



for all N e {iVo, 2^Nq, . . .}. Summing over N e {2^2^, 2^, . . .} results in 



OO 

n=l 



(3) 



E 



Are{2^2^,2^....} 



E 



7VG{2\2^,23,...} 
Af>iVo 



l(3) 



<7Vo+ E 10ct"2T^3 7v(3- 



<5ff*Tln(Ar)/4) 



< OO. 



N=No 



This completes the proof of Lemma |4.12| 

Lemma 4.13. Assume that the above setting is fulfilled. Then we have 



□ 



OO 



< oo. 



(71) 



Proof of Lemma 4-13 Define a filtration J; , Z e {0, 1, ... , Id(iV) — 1}, through 



■N 



:= o•(e"'^ keN,ve {\d{N) ~ I, Id(iV) - / + 1, . . . , ld(7V)} ) 



(72) 



for all / G {0, 1, . . . , Id(iV) — 1} and for every N G {2\ 2^, 2^, . . .}. Moreover define an J"/S(M) -measurable mapping 
Ln: n ^ {0,1, . . . , Id(iV) - 1} through Ljy := ld(V) - for every N e {2^2^23, . . .}. Rewrite these random 
times as 



Ln = ld{N) -Ln ^ ld{N) - max 



{l}u{;e{l,2,...,ld(7V)}: 3fce{l,2,...,f}: > 2iT-i }) 



i(^{\diN) - 1} U e {0, 1, . . . , Id(iV) - 1} : 3 fc e {1, 2, . . . , ^^j^}: \(('^(^)-^^-'^ 



(Id(N)-i) 

> 2 4' ' t 



for aU N e {2^, 2^, 2^, . . .}. We observe that L AT is a stopping time with respect to the filtration J; , I G 
{0,l,...,ld(iV) - 1}, for every N e {2^, 2^, 2^, . . .}. Consequently the cr-algebras for N e {2^, 2^, 2^, . . .} 
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are well-defined. By definition of (ItSI) the random variables A: E N, are independent of for every 



iV e {2\ 2^, 2^, . . .}. In particular, we obtain 

p ^(ijv-i),*; ^ ^ 



(73) 



P-almost surely for all A e B[R), fc e N and all N e {2\ 2^, 2^, . . .}. Combining ([73]) and the fact that both 
ryjv : fi — >• [0, cx)) and Ljq : — >• {1, 2, . . . , Id(iV)} are measurable with respect to J"^ then shows 



< 



jv-l).fe 



J"; 



■AT 



^7^ 



TV 



2(ijv-l) 



k=l 

N 
2(i«-i) 



< 



4(-2'""~^')2^^r-i 



N 



N 



P-almost surely for all N £ {2^22,2^,...}. Thus the inequality P[||^°'i| - a;| < e] < 4ea^'^ for aU x e M and aU 
£ £ (0, oo) yields 



< 



N 
2(i«-i) 



. 4 . 4(-2<-"-^')2^r-iiV . < ^2H"") 



(74) 



-almost surely for all N E {2^, 2^, 2"^, . . .}. Now we apply inequality (74| to obtain 



A 



< 



E 



(75) 



lw-eA'l<4H''""')2^^r-i7V 



< E 



2(-2^") 



for ah N e {2^2^, 2^, . . .}. Next we observe 



2^" > 2L21d(ff'T5 1„(Ar))J ^ 2(21d(5'T31n(W))-l) ^ 1 ((3^T 5 In(JV))^) ^ l^'^T (ln(iV))2 

on (A^^)^ for aU iV e {2\ 2^, 2^, . . .}. Inserting ([76]) into ([75]) results in 



(76) 



A 



< E 



1, ,,,^.1^2(-^-^W))^) 



< l^exp('-ln(2)-a4T(ln(7V))2') = 4 ^-i 7v(2-M2)*^Ti„w/2) 



(77) 



for ah N e {2^2^, 2^, . . .}. Combining ([77|, Lemma 4.9 and Lemma 4.10 



then shows 



^(1) 



n=l 



Ep[4t)]=EPK"n(ArJ)^n(A 

n=l 

<£p[AWn(A(^J)^n(A 

n— 1 n— 

00 

< ^ 4a-ir-i7v(^-i"(2)<.-^'"W/2) ^^pj^(2)j _^^pj^ 



^(1) 

^2" 



n=l 



+ ^p[4^„)uA« 

n=l 



n=l 



< 00. 



Af=l ri=l 

This completes the proof of Lemma |4.13| 

Lemma 4.14 (Almost sure finiteness of A'^i). Assume that the above setting is fulfilled. Then we have 



□ 



iVi < 00 



= 1. 
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Proof of Lemma \4-14\ Combining the subadditivity of the probabiHty measure P and Lemmas |49l [ilol [4T2] and 
14.131 shows 



n=l 



n=l 



< CX). 



The lemma of Borel-Cantelli (e.g., Theorem 2.7 in [23]) therefore iniphes 



(78) 



hmsup (a^I? U A^^} U A^t' U A^t' 

11, — ^OO ^ 



(79) 



Hence, we obtain 

p[7Vi < oo 



e 17 : 3 71 e {TVo, 2iiVo, 2^iVo, . . .} : Vm e {n, 2^^, 22^, . . . } : ^ uf^jA^ | 
|tj e r2; 3n e N: Vm e + 1,...}: cj ^ A^li U yl^^i U A^^i U A^t^ | 

(■ 



hminf(AWuA^^Ju4t'u4^) 



= 1. 



This completes the proof of Lemma |4.14| 



(80) 



□ 



Proof of Theorem \4.1\ Fix p £ (0,cxd) throughout this proof. Our proof of Theorem 4.1 is then divided into 
four parts. In the first part we analyze the behavior of the multilevel Monte Carlo Euler approximations on 
the events {tjn > n {Ni < N} = {lu G n-. rjNiuj) > 2(-^"('^)+i)/4r-i/4^ ^^(^) < ^| fo^, ^ g 

{2^,2^,2'^,...} (see inequality ([si])). In the second part of this proof we concentrate on the events {9n > Vn} 
n {riN < 2('^"+i)/4r-i/4| p < jsfj fQj. g {21,22,23,...} (see inequality ([91])). In the third part of this 
proof we investigate the events {2^N/4rp-i/4 ^ g^j ^ ^ ,^^| ^ < ^(Lr,+i)/4rj.-i/4^ ^ < ^| 

N G {2"'^, 2^, 2"^, . . .} (see inequality ([96|) and in the fourth part we analyze the behavior of the multilevel Monte 
Carlo Euler approximations on the events {2^n/4t-i/4 > g^j ^ ^ ^ |,^^ < 2(^"+i)/'*r-i/4} n {A^i < N} 
for A'' e {2^2^,2^,...} (see inequality (l97|). Combining all four parts (inequalities (Isl]), ([91]), ([96| and ([9?])) and 



Lemma |4 . 1 4| will then complete the proof of Theorem |4.1| as we will show below. In these four parts we will frequently 
use 

{N, <N}c {A^^^Y n {A^^^r n (Aj^V n (^j^y (8i) 

for all iVe {21,22,23,...}. 

We begin with the first part and consider the events {-qN > 2'^^'^+^'>/'^T-^^^}n{Ni < N} for N e {2^ ,2"^ ,2^ , . . . }. 
Note that Lemma [Zj] the inequalities tjn > 2^l^+^'>/'^T-^/'^ {l + 5(-*-2''")) on {tjn > 2(^"+i)/4r-i/*} n {A^i < N} 
(see ([43])) and < 2('-i)/4r-i/47v on {iVi < N} for aU fc G {1, 2, . . . , f }, ? e {1, 2, . . . , ld(7V)} (see ([42])) and the 



definition ( 38 ) of 


imply 


;=0 fe=i 


2Ljv 

> 




> 


2^«,(2'^ 

2^2^N 



i=l fc=l 



E]\^Eki 

;=1 fc=i 



2('-i).^'. 
2('-i) 



E 



1=1 



E ]VE|2/: 

i=Ljv + l fe=l 
i'2('-i) 



2(/-l) 
2('-i)' 



P Id(Af) 

- E 



(82) 



2<'-i',(#)' 



on {ryjv > 2(^"+i)/4t-i/'*} n {A^i < iV} and Lemma 



4.8 



Lemma 



4.4 



and Lemma 



4.2 



hence yield 



1=0 k=l 



y -T 

1=1 k=i 



y2('-i)' 



> 



1 

TV 



2'2 ™ \* /'r(l-*)2''« ' 



-E 



1=1 



2(1-1) 
T 



N 



(.a-D) 



p-5 



> N-^T-i -expl^l -5(1- 



-5)2^ 



- L 



Id(Af) 

E 

i=Ljv + l 

- Id(iV) 
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on {rjN > 2(^"+i)/''T-i/4} n {^i < N} for all N € {2\ 2^, 2^, . . .}. Therefore, we obtain 



=0 fc=l 



2 ,f 
^2' 



1=1 k=l 



y2('-i)' 



> 7V-ir-S . exp(| . 5(i-^)2-«^ _ ld(iV)iv5T-5 . ivl'' J - ld(iV)iv5T-5 

> T-? • exp(2 . 5(i-*)2"" - In(iV)) - T"? • iV^ 



(83) 



l+f +P-5 



on {77JV > 2(^"+i)/'*T-i/''} n {TVi < iV} and the estimate 2^" > a^VT\n{N) on {iVi < TV} (sec ^) hence shows 



H TV 5^ 

;=o k=i 



2',£' 

2/2' 



^ iV^ 

i=l fe=l 



(84) 



> inf 

i:e[ff^\/Tln(JV),oo) 



(^E . 5((i-5)^) _ in(7v)) - exp(ln(Ar) (l + f + p . 



exp § • 



55' 



T-^ > r{N) -T^S 



on {riN > 2(^"+i)/4t-i/4} n {TVi < iV} for aU N e {2^,22, 2^, . . .} where r : N ^ M is a function defined by 



r{N) := inf 



exp 



(2 . 5((i-*)-) _ ln(2A^)) - exp(ln(iV) (l + f +p • 5?) ) 



for all N eN. 

In the next step we analyze the behavior of the multilevel Monte Carlo Euler approximations on th e ev ents 
{On > Vn} n {rjN < 2(^"+i)/4r-i/4} n {Ni < N} for N e {2\ 2^,23, . . .}. To this end note that Lemma [l5| the 
inequalities 61a, > (l + 4("2*''"""))ryAr on {On > r/Ar}n{iVi < N} (see ^) and < 2('-i)/4y-i/4jy < ^} 

for ah e {1, 2, . . . , f }, / e {1, 2, . . . , Id(iV)} (see ^) and the definition ([38]) of Ljv imply 



i=0 fc=l 



2' 



E ]^Ek: 



2(1-1)^^1, 
2('-i)' 



> 



> 



iV 
1 

N 



2(-£'jv-i). 



y2(i]v-i) 



Ljv-l 



y2(i]V-i) 
2(^«-i),(l+4(-^*^""")) 

y2(i]v-i) 



Vn 



y -T 

1=1 fc=i 



W2(^iv-i) 



2/2(i-i)' 



ijv-l 

E 



Id(A') 

y -Y 



2^-1) (- 2(' -1) ^)4 

y2('-i) 



P Id(Af) 

- E 



2('-i),(|l)^ 
y2('-i) 



on {e'Ar > t^at} n {r^AT < 2(^«+i)/47^-i/4} n {A^i < iV} for aU N e {2\ 2^,23, . . .}. Lemma 
Lemma Sil] therefore show 



4.5 



Lemma 



4.4 



and 



id(JV) 

E ]vE 

!=0 fc=l 



2'.f' 
2/2' 



> 



1 

2]V 



ld(W) 2i 

E ]vE 

1=1 k=i 



y2('-i) 



2(-«-i),(l+4(-^'"""')) 

y2(-^N-i) 



y2(-^«-i) 



p 1 
2iV 



2(^"-i',7;jv 
2/2(iiv-i) 



(85) 



E ^ 



T 



iV 



p-5 



(2('-l)) 



Id(Ar) 

E 



2/\ 4 

T 



on {e^AT > 7]n} n {r/w < 2(^"+i)/4y-i/4| p < ^} ^j. N e {2\ 2^, 2^, . . .}. By definition of tjn and of Ln 
we have tjn > 2^'' I^T'^I'^ on {A^i < iV} (see (|4l])) for all iV e {2^2^,23, . . .}. Consequently we get the inequality 



T]N > 2^"/4T-i/4 (^1 + 5("*'2*''" on{7Vi < iV} (see ^) for aU TV e {2^, 2^, 2^, . . .}. Lemma 4.6 
hence yield 



and Lemma 4.5 



ld(W) „i W 

\ ^ 2 

7vEp2 

1=0 k=l 



ld(W) 2; 

E ]V E k: 

i=l fe=l 



2(1-1) 
2('-i>' 



> 



Lfi + 4H<^"-') 
2iV V 



2<^«-i',W 
"2(^iv-i) 



1 

2iV 

i^Af 



2(^-^-1). ( 2^w 



(86) 



I ^ — I N\ J - Id(iV) I 



V J" 
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on {9N>m}r\ {m < 2^^«+^'>/'^T'^/*} n {Ni < N} for all N e {2\2^ ,2^ , . . .}. Lemma [Is] and Ln < \d{N) 
therefore imply 



2', 

/=0 fc = l 



;=1 fe=i 



> 



2iV 



,2(^«-i),^„ 
^2(-f'N-i) 



(^5(..-«--))-Mm(^)VK'"""'' 



2iV V T 



exp 



id(7V)7V3r 



on {6'Ar > 77Ar} n {?7Ar < 2(^«+i)/4t-i/4| p {^v^ < TV} for all iV e {2^,22, 2^, . . .}. The inequalities 1 < Lat < Id(A^) 
and 1 + 2^^^'' > exp(2'^^^^^'') for all x £ [0, 00) hence give 



TV^ 

1=0 k=i 



2 ,f 

^2' 



y -T 

1=1 fc=i 



2/2(1-1)' 



> 



2(-C-N-i) 



(i^.exp(2(— ).p.5(^<--"))-l) 
+ ^ . exp(| . 5((i-^)2<^"-^')) - 21d(7V)7Vf T-5 • ^(^•^^^''""'^) 



(87) 



on {e^AT > m)^{m < 2'^^^+^^/^T-^/^} n {A^i < N} for an e {21,22,23,...}. This shows 



/=0 fe = l 

> ( W • exp 



1 



2' 



i=l A:=l 



P 



£.(|)(2'-"-^') 



2(1-1) ^l.fc 
2('-i) 

2('^«-i),^„ 



2ivr 

> ( w • exp 



F • exp 



Ij P2<^N-1) 
|^£.5((l-5)2<^«-i)) 



T-4 . TV 



l + f+P-5 



(2(-f'IV-2)) 



2(^~-i),^„ 
»2<^N-i) 



inf 

a:e[2(-^iV-i),oo) 



exp 



(1 . 5((i-5):^) _ ln(2Ar)) - exp(ln(iV) ( 



1 + f +p-57~ 



on {On > i]n} Ci {rjN < 2(^"+i)/4T-i/4} n {Ni < N} for aU TV e {2^, 22, 2^, . . .} and, using the estimate 2(^"-i) > 



a^VT\ii{N) on {A^i < N} (see (41)), 



E TvE 

;=o k=i 



^2' 



E TvE 

1=1 k=l 



^2(1-1)' 



»2<^iV-i) 



inf 

Ke[CT2yTln(7V),oo) 



exp(f • 5«i-^)"' - ln(2iV)) - exp(ln(Ar) (l + f + p 

on {e^AT > r;Ar}n {t]n < 2(^"+i)/4T-i/4} n {Ni < N} for ah N e {2^, 22, 2^, . . .}. It follows from 

lim ^ expf 2 . Tv(*'^i"(5/4))) = 00 

that there exists an N2 G {2^, 22, 2^, . . .} such that 

^exp(|.7v(*'^i"(^/4))) -1>0 

for all N e [N2,oo). Using this, we deduce from 



(89) 



(90) 



TV ^ r2' 

;=o fc=i 



2(1-1)^^1. 



TV /-^ |y2('-i) 
i=i fc=i 



(91) 



> 



(_^.exp(|.M^ 



'VTln(5/4))\ _ 



1 



y2<^N-i) 
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on {61 AT > tjn} n {rjN < 2'^^"+^^/^T-^/^} n {A^i < N} for all N G {N2, 2^N2,2'^N2, ...}. 

Next, we analyze the behavior of the multilevel Monte Carlo Euler approximations on the events |2^™/'*T~^/'* < 
On < VN < 2(^"+i)/'*r-i/4} n{Ni < N} for N G {2^, 2^, 2^, . . .}. Note that Lemm a [iisj and the inequality 
|^/,/c| ^ 2('-i)/4r-i/47v on {A^i < N} for aU A; e {1, 2, . . . , f }, / G {0, 1, 2, . . . , lci(7V)} (see imply 



2^2'' 



TV ^ , 

;=o k=i 

N 

> y -T 

1=1 fc=i 



;=i fe=i 



2(1-1) ^^1 
2/2(!-i)' 



i=0 fc=l 



1 

> 

- N 



tf2(iiv-i) 



> 



2iV 



»2(^iV-i) 



2(ijv-i)^ 

y2<-^«-i) ' 



2(i-N-i)^fl 
2/2(-f-jv-i) 



P X — ^ X — ^ I 2' 

/=0 fe=l 



p v^' 2 \ ^ I 2 

l — Ljsr k—1 



p 1 

^ 27V 



W2(^N-i) 



E 



2^2' 



P Id(Af) 

-E 



2'.(^)^ 
2/2' 



on {2^"/4r-i/* <0n <Vn < 2(^"+i)/4t-i/''} n {TVi < N} for all iV e {2\ 2^ 2^, . . .}. Therefore Lemma 
inequality rjN > (l +4("^*^" ^^^)0n on {6'Ar < tin} n {A^i < iV} (see (|44|) and Lemma result in 



(92) 



4.5 



the 



E ]vE 

i=0 fc=l 



2/2'' 



E ]vE 

/=! fc=l 



2/2(1-1) ' 



> 



2Ar 

Ljv— 2 

E 

(=0 



2(-«-i),(l+4(-^'"""'))^ 
2/2(ijv-i) 



P 



2<-f'N-i)_ 
2/2(iN-i)' 



1 

2iV 



y2(-f'iv-i) 



2/2' 



-E ^ X 



4.6 



(93) 



Lemma 



4.5 



on {2^"/4r-i/4 < 61a, < 77JV < 2('^"+i)/4r-i/4} n {A^i < iV} for all N e {2^,22, 2^, . . .}. Lemma 

and the estimate tjn > 2^"/'*T"i/'* (l + 5("'^'^*''""'')^ on {Ni < N} (see ^ and Q) and Lemma 4.4 hence yield 



E ]vE 

Z=0 fc=l 



2^ f ^ ' 
2/2' 



E TvE 



2/2('-i)' 



> 



AT V / 



p-5 



(2<^N-l)) 



2N 

1 

2iV 



2(^iv-i),(2£«)I(i+5-^-2(^-iv-i)^ 
^2<-^JV-i) 



- 1 
P Lm~2 



2(-f'JV-i)_ 
2/2(ijv-i)' 



E 



TV 



p-5 



- 23 ld(iV)iV3T"3 



on {2^«/4t-i/4 < On < m < 2(^"+i)/4t~i/''} n {iVi < iV} for all iV e {2\ 2^, 2^, . . .}. Therefore Lemma 
implies 



4.8 



1=0 k=l 



2/2'' 



1=1 k=i 



2'- 

2/2(!-l)' 



1 

2iV 



2<-C'iV-i),fl 
2/2(iiv-i) 



(94) 



(l + 2H""))^'^ J -I 

+ 25v J • - 21d(^)^f T-f ^(-^'^"^^'^) 

on {2iiv/47^-i/4 <9n <Vn < 2'^^"+^y'^T~^/^} n {TVi < N} for all iV e {2\ 2^, 2^, . . .}. The inequality 1 + 2'-^) > 
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exp(2'^ ^ for all x € [0,oo) hence shows 



1=0 k=l 



2/2'' 



y -T 

1=1 k=l 



> 



,(|.5((i-^)2''""'))-T-?.7V 



y2(ijv-i)' 



(95) 



27VTt 



F ■ exp 



l+f +P-5 



on {2^N/4rp~i/4 < 6»jv < ryw < 2('f'"+i)/4T-i/4| ^ < ^} foj. iV e {21,2^,23, . . .}. Consequently 



;=o k=i 



1 



y2'' 



1=1 k=i 



^(^•-P(f-(l) 



_ -, I I 2<^"-i>, 
-l- ) |2/2(i-N-i) 



inf 

a;e[2(-^iV-i),oo) 



exp(| • 5«i-*)") - ln(2iV)) - exp(ln(iV) (l + | +p • 5^) ) 



on {2^"/4t-i/4 < 6'Ar < r^AT < 2(^"+i)/4T-i/4} n {A^i < TV} for all N e {2\ 2^, 2^, . . .}. The estimate 2(^"-i) > 



CT2VTln(7V) on {Ni < N} (see (41)) therefore implies 



^(2]V-°"P 



2' 



2('-i).5'.' 



E_ \ ^ 2 

Z=l fc=l 



2 ' U; 



- 1 



y2(^N-i) ' 



inf 

i:e[CT^\/Tln(7V),oo) 



exp 



(^E . 5((i-5)^) _ in(2iV)) - exp(ln(Ar) ( 



l + f +p-5^ 



on {2^"/4r-i/4 < 6'jv < W < 2('^"+i)/4r-i/4} p < ^j. foj. all iV e {2\ 2^, 2^ . . .}. Finally, we obtain 



/=0 fe=l 



2' 



- E ]^ E k: 

Z=l fe=l 



2('-i) 



(96) 



> 



2N 



■ exp 



/Tln(5/4)) 



2(-^«-i), 
-l- ) |j/2(-f'N-i) 



+ r(iV) • T-^ > r{N) ■ T^i 



on {2^"/4r-i/4 < ^jv < ?7Ar < 2('f'"+i)/4T-i/4} n {Ni < N} for all N € {N2,2^N2, 2^N2, . . .}. 

Finally, we analyze the behavior of the multilevel Monte Carlo Euler approximations on the events {9^ < 
2Liv/4y-i/4| p <rjN < 2^LM+i)/4rp-i/4^ P < ^| f^j. ^ ^ {21,22,23,...}. Note that Lemma 



4.5 



and the 



inequality j^'^'^l < 2('-i)/'^T-i/47V on {TVi < N} for aU fc e {1, 2, . . . , f }, / e {0, 1, . . . ,ld(7V)} (see ^) imply 



I 2',f''''' 
^2' 



E ]vE , 

N 

- E ]V E |2^2('-i) 
i = l fc=l 



E TV E 

; = 1 A;=l 



2/2('-i)' 



1=0 k=l 



^2^1 .k 
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> 
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1 

N 



y2(-f-]v-i) 



Ln~2 2^ 



1=0 k=l 



2',C' 



]^EP2' 

/=Ljv-l fc=l 



2(-f.iv-i) 



P ijv-2 



E 

(=0 



^2' 



ld{N) 

E 



2^(^)^ 
i/2' 



on {6'iv < 2^"/4T-i/4} n {On < W < 2(^"+i'/4t-i/4} n {TVi < N} for all iV e {2^, 2^, 2^, . . .} and, applying 
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Lemma 14.81 Lemma 14. 51 Lemma 14.41 and Lemma 14.21 



1=0 k=l 



2^ 7 ^ 

^2' 



2(1-1) ^i,k 



/=1 A: = l 



> 



TV 



2ijv\4 ^5(1-5)2(^1V-1) 



P Liv-2 /2l\i 
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i=0 



T 



\d{N) 

E 

=L„-1 



2 • 2 



/ \ 4 



> TV^^T^i -exp 



on {0N < 2^"/4r-i/4} n {6'jv < Vn < 2(^"+i)/4r-i/4} n {iVi < N} for all iV g {2\22,23,...}. Therefore, we 
obtain 



y -T 

1=0 k=l 



y2'' 



^ N ^ 

1=1 k=i 



> iV^^T^i • exp 



-5)2(^N-i) 



p-5 



> T-? • exp(f • 5(i-^)2<""-^' - In(iV)) - T"? • TV 



ld(iV)iV3T"3 • N 



on < 2^"/'*T ^Z'*} n I^AT < ?7Ar < 2(^™+^)/^r ^/^} H {iVi < N} and hence, using on 



ivEp2'' 

/=0 fc=l 



> inf 



exp 



E ]^ E 

i=i k=i 

(^E . 5((i-5)^) _ in(2TV)) ~ exp(ln(7V) (l + | +p • 5?) ) 



T"4 ^ r(/V) • 



(97) 



on {On < 2^«/4r-i/4| n < t^at < 2(^"+i)/4t-i/4} n {A^i < iV} for all N e {2\2^,2^, . . .}. 
Combining (|84|, ml, fctl) and 1971) then shows 



E ]^Ek: 

/=0 k=l 



2' 



p ld(W) ^ 

- E ]^ E k: 

i=l fc=l 



2(i-l) ^^fc 
2('-i)' 



> r(7V) -T-i 



(98) 



on {A^i < iV} for aU N e {iVa, 2iiV2, 22A?^2, • ■ • }• Equation (37) and inequality (98 1 imply 



1 ld{N) F / X 

^2^|n ' (W)| + 2^ ]^ 2^(1^2' HI -F2('-i) HI ) 
*:=! 1 = 1 k=l ^ ^ 



^ iV ^ 

1=0 k=\ 



I 2'.5'''=(w) 



F -Y 

^ N ^ 

1=1 k=l 



I 2<'-i',5''''-(a.) 
P2('-i) 



> r(iv) • r-3 



(99) 



for all N e {iVi(a;),2i • A^i(tj),22 • iVi(w),...} n [iV2,oo) and all w e {TVi < oo}. The fact liniAr^oo r-(iV) 
therefore shows 



lim 

\d(N)m 



1 N Id(JV) iT / s 

]VEI^1 HI +Z^]^Z^( 1^2' HI -|i"2a-i, (w)| 1 

k=l 1 = 1 k=l ^ ^ 



OO 



(100) 



for all u) e {A'^i < oo}. Hence, Lemma 4.14 finally yields 



lim 

ld(Af)eN 



1 ^ 

-Fin' 



1,0, fc IP 



fc=l 



;=i k=i 



2',Lk\P _ |y2<'-i',/,/e| 



'2(1-1) 



(101) 



-almost surely. This completes the proof of Theorem |4.1 



□ 
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5 Divergence of the multilevel Monte Carlo Euler method 



Motivated by Figure |4] below and by the divergence result of the multilevel Monte Carlo Euler method in Section |4j 
we conjecture in this section that the multilevel Monte Carlo Euler method diverges with probability one whenever 



one of the coefficients of the SDE grows superlinearly (see Conjecture 5.11. Whereas divergence with probability 
one seems to be quite difficult to establish, strong divergence is a rather immediate consequence of the divergence 



of the Euler method in Theorem 2J_ above. We derive this strong divergence in Corollary |5.2| below. For practical 
simulations the much more important question is, however, consistency and inconsistency respectively; see, e.g., 
Nikulin [53], Cramer [2 , Appendix A.l in Glassermann [TO] and also Theorem 4.1 above and Conjecture 5.1 below. 



Throughout this section assume that the following setting is fulfilled. Let T S (0, oo), let (fi, P) be a probability 
space with a normal filtration {J^t)te[o,T], let W^''' : [0, T] x — M, / G No, A: G N, be a family of independent one- 
dimensional standard (J-t)tg[o^T]-Brownian motions and let — > M, ^ e No, /c G N, be a family of independent 
identically distributed Jo/'B(R)-nieasurable mappings with E[|^°'^|p] < oo for aU p G [1, oo). Moreover, let fi,a: M ^ 
R be two continuous mappings such that there exists a predictable stochastic process X : [0, T] x f2 ^ M which satisfies 

|/i(Xs)| + |cr(Xs)|^ ds < oo P-almost surely and 



Xt^S.°^'+ / / a{X,)dW^'^ (102) 

Jo Jo 

P-almost surely for all t G [0,T]. The drift coefficient /i is the infinitesimal mean of the process X and the diffusion 
coefficient a is the infinitesimal standard deviation of the process X. We then define a family of Euler approximations 



YN,l,k. Q 



n G {0, 1, . . . , iV}, TV G N, Z G No, /c G N, by 



N,Lk 



^'■'^ and 



-'n+1 ■— 



N.Lk\ 



Lk 



w: 



Lk 



for aU 71 G {0, 1, . . . , iV - 1}, G N, / G No and all fc G N. 



Conjecture 5.1 (Divergence with probability one of the multilevel Monte Carlo Euler method). Assume that the 
above setting is fulfilled and let a,c € (l,oo) be real numbers such that ^J— < |/i(a;)| -I- |cr(a;)| < c\x\'^ for all 
X £ M. with \x\ > c. Moreover, assume that P[cr(^'''^) 7^ O] > or that there exists a real number (3 G (l,oo) 

such that P[|^°'^| > a;] > ^("^'') for all x G [l,oo). Moreover, let f : R ^ M. be B{R)/B{R) -measurable with 
-\x\^ — c < f{x) < c(l + jxl"^) for all x eM.. Then we conjecture 



lim 

ld(Af)GN 



N 



E 



N 



1.0, k\ 



k=l 




k=l 



(103) 



V-almost surely. 

To support this conjecture, we ran simulations for the stochastic Ginzburg-Landau equation given by the solution 
(-'^t)te[o,i] of 

dXt = (2Ai - Xl) dt + 2Xt dWt, Xo = 1 (104) 
for all t G [0, 1]. Its solution is known explicitly (e.g. Section 4.4 in [55]) and is given by 

exp {2Wt) 



Xt - 



1 + 2 exp {AWs) ds 



(105) 



for t G [0,1]. We used this explicit solution to compute E[(Ai)^] « 0.8114. Figure |4] shows four sample paths of 
the approximation error of the multilevel Monte Carlo Euler method for the Ginzburg-Landau equation (104). Only 



finite values of the sample paths are plotted. The next corollary is an immediate consequence of Theorem |2.1| above. 



Corollary 5.2 (Strong divergence of the multilevel Monte Carlo Euler method). Assume that the above setting is 
fulfilled and let a,c £ (l,oo) be real numbers such that < |/^(a;)| + |cr(x)| < c\x\'^' for all x £ R with \x\ > c. 
Moreover, assume that P^a{£,^'^) 7^ O] > or that there exists a real number (3 G (l,oo) such that P[\£_^'^\ > a;] > 

for all X G [l,oo). Additionally, let f : R ^ R be B{R)/B{R) -measurable with ^\x\- -c< f{x) < c(l -I- jxl'^) 
for all X G M. Then we obtain 



lim 

ld(7V)eN 



E 



N 



E 



fe=l 



for all p G [1, 00). 




>(i-i) 

2(i-i) 



,Lk^ 



fe=l 



(106) 
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Proof of Corollary \5.^ First of all, note that the assumption E[|^°'^|p] < oo for all p e [l,oo), the continuity of 
/I, a: M ^ M, the inequality |/i(a;)| + < c \x\'^ for all a; e M with > c and the estimate \f{x)\ < c(l + for 

all a; e M imply E[|/(Y^'"'^)|] < oo for aU N eN. Therefore, we obtain 



E 



N 



N 



fc=i 



1=1 \k=i 



= E 



N 



for all N e {2\ 2^, 2^, . . .}. The estimate f{x) > l\x\^ - c for all x e M and Theorem 2.1 hence give 



lim E 

ld(Af)eN 



N 



^E/(n^'°"')+ E ^7 E/(^^'"v/(i^.t; 



k=l 



i=i \fe=i 



) J 



> - hm E 

C VAf-s-oo 



r, 



iV,0,l 



In the case E[|/(X7-)|] < oo, the triangle inequality and Jensen's inequality then yield 



hm 

Af-i-oo 

id(Ar)eN 
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k=l 



> hm E 

Ar->-oo 

ld(Af)eN 
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;=i \k=i 



TV 



A:=l 



1=1 \k=l 



E 



/(^t)| 



oo 



(107) 



forallpS [l,oo). This shows ( 106 ) in the case E [|/(Xt)|] < oo. In the case E = oo, the estimate /(.t) > — 

for all X S M shows E[/(Xt)] = oo and this implies (106) in the case E[|/(Xt)|] = oo. The proof of Corollary 
is thus completed. 



5.2 



□ 



6 Convergence of the multilevel Monte Carlo tamed Euler method 

In this section we combine the multilevel Monte Carlo method with a tamed Euler method. We aim at path-dependent 
payoff functions. Therefore, we consider piecewise linear time interpolations of the numerical approximations, which 



have continuous sample paths and which are implementable. Thereorem 6.1 shows that these piecewise linear 
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interpolations of the tamed Euler approximations converge in the strong sense with the optimal convergence order 



according to Miiller-Gronbach's lower bound in the Lipschitz case in Theorem 6.2 then establishes almost 

sure and strong convergence of the multilevel Monte Carlo method combined with the tamed Euler method. The 
payoff function is allowed to depend on the whole path. We assume the payoff function only to be locally Lipschitz 
continuous and the local Lipschitz constant to grow at most polynomially. 

Throughout this section assume that the following setting is fulfilled. Let T G (0, cx)), let (17, P) be a probability 
space with a normal filtration (~^t)te[o,T], let d, m g N, let W'-'^: [0,r] x O M™, ; e Nq, fc e N, be a family of 
independent standard (J"t)tg [o.T]-Brownian motions and let : Q, — > W^, I e Nq, fc e N, be a family of independent 
identically distributed Jx)/S(M'')-measurable mappings with E[||^°'i||^^] < CO for all p G [1, oo). Here and below we 
use the Euclidean norm |ja;||]g„ ■— + x\ + ■ ■ ■ + x^-^ for all x = (xi, a;2, . . . , x„) e M" and all n e N. Moreover, let 
/i: M'^ — ^ M'' be a continuously differentiable and globally one-sided Lipschitz continuous function whose derivative 
grows at most polynomially and let a = (o'ij)ig{i.2,....(i},je{i,2,...,m} '■ j^dxm ^ globally Lipschitz continuous 

function. More formally, suppose that there exists a real number c € [0, oo) such that (x — y, ijl{x) — ljL{y))^d < 
c\\x - 2/||Rd, \\^l'{x)\\nJ^.d) < c(l + ||a;||Rd) and \\a{x) - cr(2/)||L(R'»,R'') < 4^ - vWr" for all x,y e R'^. Here and below 



we use := {J2Li \^^?)^' 
consider the SDE 



and {x,y)g^i := ' Vi for all x = (xi,X2, ■ . . ,Xd),y = {yi, ■■■,yd) & R''- Then 



dXt = ^liXt) dt + (j{Xt) dWt 



0,1 



^0 



0,1 



(108) 



for t e [0, T]. Under the assumptions above, the SDE ( |108 ) is known to have a unique solution. More formally, there 
exists an up to indistinguishability unique adapted stochastic process X : [0, T] x 17 — > M'* with continuous sample 
paths fulfilling 



Xt^C' + / fi{Xs)ds+ / aiX,)dW, 



^0,1 



(109) 



P-almost surely for all t S [0,2^] (see, e.g.. Theorem 2.4.1 in Mao [IH])- The drift coefficient ^ is the infinitesimal 
mean of the process X and the diffusion coefficient a is the infinitesimal standard deviation of the process X. In 
the next step we define a family of tamed Euler approximations Y^'^'^ : J7 W^, n E {0, 1, . . . , N}, N E N, I E Nq, 



fc e N, by Fo 



^'"'^ and 



N,l,k yAr,;,fe 



N.Lk\ 



T 

N 



1 



fT^N,l.k\ 
KYn ' ) 



N,Lk\ 



(n + l)T 



w: 



Lk 



(110) 



for all n E {0, 1, . . . , — 1}, N E N, I E Nq and all k E N. In order to formulate our convergence theorem for the 
multilevel Monte Carlo tamed Euler approximations, we now introduce piecewise continuous time interpolations of 
the time discrete numerical approximations ( |110[ ). More formally, let y^^'-*^ : [0, T] x O ^ M'', G N, Z e No, k eN, 

(Ill) 



be a family of stochastic processes with continuous sample paths defined by 



T 

N 



(^N,l.k 



T 



n Y 



N,l.k 
n+1 



1 



J n 



nT {n+iyr^ ^ c /n 1 - 1}, iV e N, / e No and ah fc e N. 

The following corollary is a direct consequence of Hutzenthaler, Jentzen and Kloeden 



for alHe [^,^^],ne{0,l 



^ and Miiller-Gron- 

bach [33] (see also Ritter [3S]). It asserts that the piecewise linear approximations F^, N E N, converge in the 
strong sense to the exact solution. The convergence order is | except for a logarithmic term. 

Corollary 6.1 (Strong convergence of the tamed Euler method). Assume that the above setting is fulfilled. Then 
there exists a family Rp E [0, oo), p E [1, oo), of real numbers such that 





sup 


Xt - 


p 






tG[0,T] 







Vn 



(112) 



for all N eN and all p E [1, oo). 



The convergence rate N 2 (l-|-ld(iV)) 2 for N eN obtained in ( 112 ) is sharp according to Miiller-Gronbach's lower 
bound established in Theorem 3 in [33] in the case of globally Lipschitz continuous coefficients (see also Hofmann, 
Miiller-Gronbach and Ritter [HI). 



Proof of Corollary\6^ Let Y^ : [0, T] x f7 ^ W\ iV e N, be stochastic processes defined by 



Y, 



N 



Y, 



1 



\lJ'{Yn ■ ) 



N ) 



a{Y:^''^'){w^-'-W'j) 
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for all t e ^^^7^], n e {0, 1, . . . , iV - 1} and all N e N. Theorem 1.1 in then shows the existence of a 

family Rp £ [0, oo), p e [l,oo), of real numbers such that || sup^gjo.T] W-^t ~ ^t^lk'' ||lp(o r) — ^ ^ ^ 

all p € [IjCxd). The triangle inequality hence yields 



sup 

tG[0,T] 



LP(0;I 



sup 

tG[0,T] 



(113) 



LP(0;R) 



for all e N and all p G [1, oo). Moreover, we have 



_ yNfi.l 



IN 



0,1 



"TV « 



(114) 



<||a(r„~'°^i)| 



tN 



TT/O.l tiaO,1 / ''^^ \ / TirO.l TT/0.1 



for aU < e ^""^^^^ ], n e {0, 1, . . . , TV - 1} and all iV e N. Combining ([Tl3|, (ITu]) and Holder's inequality then 



gives 



sup 

te[o,T] 



N,0,1 



< 



Rn 



N 



max sup 

rie{Oa,...,Ar-l} (,i+l)T l 



< —j^ + A/ T7 C • sup 



max 1 1 Y, 



M,0,1 II 



lk(o)l! 



L2p(ri;B 



L(R™,R'') 



L2p(a;] 



L2p(a:R) 

max sup |/3r-t-/3i"| 
xe{i,2,...,Jv}tg[o,i] 



(115) 



L2p(n;B 



for all iV e N and all p e [l,oo) where /3" : [0, 1] x 17 — > M, rt e N, is a sequence of independent one-dimensional 
standard Brownian motions. Moreover, Theorem 1.1 in [21 , in particular, implies 



sup 

Men 



max 1 1 
ie{o,i,...,J\/} 



M,o,i II 



< OO 



(116) 



LP(0;I 



for all p E [l,oo). Additionally, Corollary 2 in Miiller-Gronbach (see also Ritter [3S]) shows 



sup I (1 + ld(Ar))"3 



max sup \(3l'-t-p^\ 
7ie{os,...,N} t(z[a,i] 



< oo 



(117) 



LP(a;I 



for all p e [l,oo). Combining (115), (116) and (117) finally completes the proof of Corollary 6.1 



□ 



Proposition 6.2 (Strong consistency, converence with probability one and strong convergence of the multilevel Monte 
Carlo tamed Euler method). Assume that the above setting is fulfilled, let c £ [0,oo) and let /: C([0,r],M'^) — > M 6e 
a function from the space of continuous functions C([0, T], M'') into the real numbers M. satisfying 

ll/(") — /(^)llc([0,T],R<') < C ^1 + ||w||(^QQ J.] Rd) + ll''^llc([0,T],R'')) 11^ ~ ^llc([0,T],R'') (H^) 

for all v,w £ C([0,r],M'*). Then there exists a family Cp G [0,oo), p G [l,oo), of real numbers such that 
/ 
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N 



E 



k=l 



=1 \fc=i 



< a . il±}^ (119) 



N 



for all N € {2^,2^,2^,...} and all p G [l;Oo). In particular, there are finite J- /B{[0,oo)) -measurable mappings 
Ce'-ft^ [0, oo), e e (0, ^), such that 



E 



/wl-^E/(^"°'')- E I E/(i'^'''''^)-/(>^''"''") 



fe=l 



;=1 \fc=l 



< 



(120) 



/or all N E N and all e E {0,h) V- almost surely. 
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The convergence rate TV^s (1 + ld(7V))2 for iV e N obtained in mm is the same as in Remark 8 in Creutzie, 
Dereich, Miiller-Gronbach and Ritter [3,. For numerical approximation resuhs for SDEs with globally Lipschitz 
continuous coefficients but under less restrictive smoothness assumption on the payofi' function, the reader is referred 
to Giles, Higham and Mao 9J and Dorsek and Teichmann [5]. Moreover, numerical approximation results for SDEs 
with non-globally Lipschitz continuous and at most linearly growing coefficients can be found in Yan [42 , for instance. 



Proof of Proposition 6.2 The triangle inequality gives 
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;=1 \fe=i 



LP(n;l 
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f{x) -E 
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k=l 
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fc=l 



LP(n;R) 



for all N e {2^,2^,2^,...} and all p £ [l,oo) and the Burkholder-Davis-Gundy inequality in Theorem 6.3.10 in 
Stroock [36 shows the existence of real numbers Kp G [0, oo), p e [1, oo), such that 
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E 



'^0-1) -/(r2*Al) + j(;p2^'-'Al) 



LP(n;R) 



for all TV S {2^,2^,2^, . . .} and all p E [l,oo). In the next step estimate (1181, Holder's inequality and the triangle 
inequality show 
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(l + ll^llL2c(o.c([o,T],R<i)) + ll^^'"'"^llL2<:(a;C([0,T],R'i)) 
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/(y2^o.l)_;(y2-^^o4) 



and Corollary 6.1 and again estimate (118) hence give 
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< 2ci?2 1 + sup y 

V Men 
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v/l + ld(TV) _^2i^^ 11^^-1,0,1^ 
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^ II- - ' iiLP(n;R) 

y2',0,l _ y2<'"i),0,l 



L2p(a;C([0,T],R<i)) 



for all TV e {2^, 2^, 2^^, . . .} and all p G [l,oo). The triangle inequality, again Corollary 6.1 and the estimate 
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l/(^^)llc([o,T].K^) < (2c+ |l/(0)||c([o,T].K-))(l + IMctlW^)) for all V e Ci{0,T],R^) then yield 
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N Id(JV) . 
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llLP(<=+i)(n;C([0,T],R'')) 



cKpR2p I 1 



sup y 

AfeN 



Af,0,l I 



\L'^P'=(n:C{[0,T]M'')} 



1 

"^^2(HvrTid(20 



and finally 



E 



fix) 



fc=i 




k=l 



LP{fl;TS 



< 2cR2 1 



sup \\Y 

MeN 



M,0,1 1 



\L^'={n;C{[0,T],R'^)) 



(l + ld(7V))^ 



2ifp(2c+||/(0)||c([o 
12cKpR2p ( 1 



N 



\{c+l) 

lLp(<^+i)(n;C([0,T],R<i)) 



(l + ld(jV))^ 



sup F 

MeN 



llL2p<=(n;C([0,T],R<i)) ) /]y 



for all N e {2^,2^,2^,...} and all p € [l,oo). This shows (119 1. Inequahty (120) then immediately follows from 

□ 



Lemma 2.1 in Kloeden und Neuenkirch [25j. This completes the proof of Proposition 6.2 



It is well-known that the multilevel Monte Carlo method combined with the (fully) implicit Euler method converges 
too. The following simulation indicates that this multilevel Monte Carlo implicit Euler method is considerably 
slower than the multilevel Monte Carlo tamed Euler method. We choose a multi-dimensional Langevin equation as 
example. More precisely, we consider the motion of a Brownian particle of unit mass in the d-dimensional potential 
with d = 10. The corresponding force on the particle is then x — |la;||^ • x for a; G M''. More 



formally, let T = 1, m = li = 10, ^ = (0, 0, . . . , 0), fi{x) = a; — ||a;|| ■ x for all x e 
matrix for all x G M''. Thus the SDE ( |108[ ) reduces to the Langevin equation 



and let a{x) = I be the identity 



{Xt - \\Xt\\ 



Xt]dt + dW 



0,1 



for t e [0, 1] . Then the implicit Euler scheme for the SDE (1211 is given by mappings : fl 
N eN, satisfying = ^ and 

T / - : 2 

N 



Y, 



N 



■Y, 



N 



w 



0,1 



(121) 

£^nG {0,1,...,7V}, 
(122) 



for all n E {0, 1, . . . , A^— 1} and all N E N. Note that we used the Matlab function fsolve(. . .) in our implementation 
of the implicit Euler scheme ( 122 1. Figure[5]displays the root mean square approximation error of the multilevel Monte 

Carlo implicit Euler method for the uniform second moment E supjgjg ||Xt||^ of the exact solution of (121) as 

function of the runtime when N G {2^, 2^, . . . , 2^^}. In addition Figure [s] shows the root mean square approximation 



suPtg[o,i] \\Xt 



of 



error of the multilevel Monte Carlo tamed Euler method for the uniform second moment E 

the exact solution of (121) as function of the runtime when N E {2^, 2^, . . . , 2^^}. We see that both numerical 
approximations of the SDE (121) apparantly converge with rate close to ^. Moreover the multilevel Monte Carlo 



implicit Euler method was considerably slower than the multilevel Monte Carlo tamed Euler method. This is 
presumably due to the additional computational effort which is required to determine the zero of a nonlinear equation 



in each time step of the implicit Euler method (122). More results on implicit numerical methods for SDEs can be 
found in [H [13 im IMl 133 m Hg , for instance. 
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